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OPTIMAL THREE-DIMENSIONAL REUSABLE TUG TRAJECTORIES FOR PLANETARY 
MISSIONS INCLUDING CORRECTION FOR NODAL PRECESSION 

by Janos Borsody 
Lewis Research Center 

SUMMARY 

Maximum payload capability for a reusable Tug is calculated by using the Maximum 
Principle to derive the equations governing the problem. The analysis includes a cor- 
rection for precession of the Space Shuttle (SS) orbit. The two-point boundary value 
problem is solved by using finite difference partial derivatives and a Newton -Raphson 
iteration technique. 

Data are presented for a given planetary mission and cryogenic Tug and include, in 
addition to payload capability, several trajectory parameters that are necessary to de- 
fine the flight profile and the outgoing asymptote. 

Inclination of the initial SS orbit is varied from 28. 5° to 90° to show the effect of 
inclination on Tug performance. Other parameters varied are Tug thrust, total trip 
time, and payload -separation coast time. Total trip time is defined as the time between 
Tug departure from and return to the SS orbit. Payload -separation coast time is the 
time from the Tug outbound burn cutoff to retroburn start. 


INTRODUCTION 

At the present time the National Aeronautics and Space Administration (NASA) is 
developing a reusable Space Shuttle (SS). As it becomes operational, many unmanned 
missions will be launched with the SS. It will carry and directly deploy payloads to low 
Earth orbits. For high -Earth -orbits and planetary missions, however, an upper stage 
(carried in the SS cargo bay along with the payload) is required. One option available for 
performing these higher energy missions is to use an expendable stage in conjunction 
with the SS. To analyze the performance and launch availability of this configuration for 
planetary missions is relatively simple using methods already developed for expendable 
launch vehicles. 



However, past NASA studies indicate that in order to reduce space transportation 
and payload costs, the SS upper stage should be reusable. This new, reusable stage is 
generally called a Tug. Performance analysis of the SS and Tug transportation system 
is far more complicated than when the expendable configuration is used, since the Tug 
must be returned to rendezvous with the SS for recovery and reuse. 

This report presents a method for determining the maximum performance of a re- 
usable SS upper stage for planetary missions and demonstrates this method by analyzing 
a particular stage for several sample cases. The mathematical formulation of the prob- 
lem is three dimensional, and a two -point boundary value problem is derived by using 
the Maximum Principle. 

Methods that have been used to determine reusable Tug performance for planetary 
missions are largely limited to ideal impulsive calculations and planar integrated tra- 
jectory simulations (refs. 1 and 2). These methods do not allow for a complete definition 
of planetary mission requirements, nor do they account for nodal precession. 

The analysis presented in this report relaxes the planar constraint and allows out- 
of -plane thrusting. It also includes a complete definition of a planetary mission in terms 
of the declination and right ascension of the outgoing asymptote and the mission energy. 
The right ascension of the outgoing asymptote is not constrained in the analysis. How- 
ever, with the information provided, the necessary right ascension can be satisfied by 
proper selection of SS launch time. 

The flight profiles presented are designed to correct for nodal precession of the SS 
orbit during the mission. Nodal precession (nodal shift) is caused by the Earth’s oblate - 
ness, which induces a retrograde motion of orbital planes about the polar axis. This 
rate of rotation is about 0. 3 degree per hour for low Earth orbits. Since the time be- 
tween Tug departure from and return to the SS orbit may be as long as 2 days, nodal 
shift becomes a major concern and a trajectory consideration. 

The maximum payload capability of a representative Tug configuration is presented 
in this report for a range of declinations and a range of departure -orbit inclinations. 

The mission energy selected corresponds to a typical inner -planet delivery mission. Ef- 
fects of total trip time, Tug initial thrust -to -weight ratio, and payload -separation coast 
time are also discussed. 


ANALYSIS 

Assumptions and Simplifications 

A simplified trajectory and flight profile for the mission to be analyzed in this report 
is presented in sketch (a). It consists of an outbound burn, a retroburn, and apogee and 
perigee burns. At outbound burn completion, three conditions that define the outgoing 
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asymptote of the escape hyperbola for a planetary mission must be satisfied. They are 
the declination and right ascension of the outgoing asymptote and the mission energy. 

The outbound burn is terminated when the desired mission energy is reached. At this 
point the engine is shut down and payload separation is initiated. 

During the fixed coast phase that follows engine shutdown, the payload is separated 
and the Tug is turned around for proper alinement to perform a retroburn. Since during 
retroburn the Tug engine is pointed toward the spacecraft, there is concern about ex- 
haust impingement on the spacecraft. In order to avoid this problem, a separation ve- 
locity must be imparted to the spacecraft, and the coast phase must be long enough to 
allow for a safe separation distance between the Tug and the spacecraft. At the end of 
the payload -separation coast phase (fixed time from outbound -burn cutoff to retroburn 
start), a retroburn maneuver is performed that reduces hyperbolic excess velocity to 
subescape velocity and puts the Tug into an elliptic orbit. The length of this burn is pri- 
marily dependent on the given Tug trip time. Trip time is defined as the elapsed time 
from Tug outbound -burn ignition to return to the SS orbit. At the end of retroburn, the 
Tug is constrained to be in the same orbital plane that the SS will be in at rendezvous. 

Retroburn is followed by a coast phase to the apogee of this intermediate Tug orbit, 
where a small delta velocity is added to adjust perigee altitude to correspond to the per- 
igee altitude of the SS orbit. This velocity change is followed by a coast phase to perigee, 
where the apogee of this intermediate orbit is reduced to correspond with the apogee 
altitude of the SS orbit. 
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Since the SS remains in a low Earth orbit while the Tug delivers the spacecraft and 
returns to rendezvous with the SS, there is a substantial nodal shift and corresponding 
plane change between the SS orbit at Tug departure and the final SS rendezvous orbit. 
Nodal shift is caused by the Earth's oblateness. The added attraction of the equatorial 
bulge introduces a force normal to and toward the equatorial plane. The resultant accel- 
eration will cause the SS to cross the equator before the point at which it would cross if 
a spherical Earth were assumed. This motion of the line of nodes is called nodal pre- 
cession or nodal shift. Nodal precession is discussed in more detail in the following 
section. Since the SS does not have the performance capability to substantially change 
the line of nodes, the Tug must perform out -of -plane maneuvers to correct for nodal 
precession. In this report it is assumed that all nodal correction is done during the out- 
bound burn and retroburn. Alternative approaches would be to introduce an additional 
nodal correction burn or to combine some nodal correction with the apogee and perigee 
burns. These alternative approaches were not evaluated, and it is not clear if they offer 
any theoretical improvements over the approach used in this report. Based on the re- 
sults obtained, the added nodal correction burn, even if theoretically optimum, may not 
be justified, since the payload loss due to nodal correction is quite small (section 
RESULTS AND DISCUSSION). While small payload gains might be theoretically realized 
by the addition of a nodal correction burn, the losses associated with an additional start- 
up and shutdown could nullify these payload gains. In deriving the equations of motion 
governing the problem, the following assumptions were made: 

(1) A spherical nonrotating Earth model is assumed. This implies that there is no 
nodal shift or apsidal precession of the intermediate Tug orbits. The Tug will correct 
for SS orbital precession only. This gives somewhat conservative Tug performance, 
since if the Tug orbits were allowed to precess, total nodal correction required would be 
slightly reduced and consequently payload capability would increase. For the mission 
analyzed herein, the nodal shift of intermediate Tug orbits is less than 2 percent of the 
nodal shift of the SS orbit. Although a spherical Earth model is assumed in order to 
simplify the equations, a nonspherical Earth model could easily be incorporated without 
a major change in the analysis. 

(2) The return leg of the mission (apogee and perigee maneuvers) is done impul- 
sively. This assumption neglects gravitational and thrust -to -weight losses encountered 
during these burns. Apogee burn is very short and these losses would be very small 
and have a negligible effect on payload. Perigee burn is longer; but by this time most 
of the propellant is expended and consequently thrust-to -weight ratio is high. Thus, the 
burn is relatively short and the losses should be very small. These burns could be in- 
tegrated, but this would significantly increase the iteration size in solving the two -point 
boundary value problem and hence greatly increase the difficulty of obtaining solutions. 

(3) The cryogenic Tug used herein has a constant thrust and specific impulse, al- 
though the analysis is formulated to accommodate a variable -thrust Tug. 
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(4) The SS is injected into a circular orbit with a given inclination and altitude. The 
Tug can depart from this orbit at any point that is determined to give maximum payload. 
A circular departure orbit is selected to simplify Tug and SS rendezvous. This is dis- 
cussed in greater detail later. 

(5) The right ascension of the outgoing asymptote is not constrained in the trajectory 
optimization. A "pseudo right ascension" is defined as the longitude of the projection of 
the outgoing asymptote in the equatorial plane measured from the ascending node of the 
SS orbit. The true right ascension, then, is the sum of the pseudo right ascension and 
the longitude of the ascending node of the SS orbit as measured from the vernal equinox. 
This longitude may be varied by varying SS launch time. Hence, for a given pseudo right 
ascension, any desired right ascension of the outgoing asymptote may be obtained by se- 
lecting the proper SS launch time. Therefore, in this report an optimum pseudo right 
ascension is determined for each declination of the outgoing asymptote, and the right 
ascension requirements are assumed to be satisfied by selection of SS launch time. 

If any unaccounted-for delay is encountered after the Tug is deployed in orbit, 
achieving the correct right ascension and declination becomes an important constraint in 
determining a Tug "launch window. " This problem is not analyzed in this report. 

(6) The analysis is done for a two -body problem, and all trajectory perturbations due 
to the Sun, Moon, and other planets are neglected. 


Nodal and Apsidal Precession 

Nodal precession has been analyzed and discussed in reference 3. An analytical ex- 
pression of secular nodal shift using first-order expansion in the dominant oblateness 
parameter J is given by the following equation: 



(Symbols and notation used are defined in appendix A. ) Nodal precession is a function of 
orbital inclination, eccentricity, and semilatus rectum. It also depends on the length of 
time spent in orbit T D (in this case total trip time), the equatorial radius and gravita- 
tional constant of the Earth, and the coefficient of the second harmonic J in the gravita- 
tional potential function. Equation (1) is used in this analysis to compute nodal preces- 
sion for a given mission. 

The line of apsides also advances on the celestial sphere as noted in reference 1. 

An analytical expression for this motion is given in terms of the change in argument of 
perigee as follows: 
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(2) 



This apsidal motion becomes very important if elliptic SS orbits are used, since proper 
apsidal orientation must be achieved between the Tug and SS orbits. For circular orbits, 
the apsidal precession translates into an effective reduction in the period of the SS orbit. 
However, this can easily be taken care of with minor payload loss or gain by varying the 
total trip time a few minutes. Although, theoretically, proper phasing of Tug and SS can 
be achieved by varying the total trip time, the operational aspects of this approach on 
Tug guidance hardware, software, and mission design requirements must be carefully 
evaluated. However, this is outside the scope of this report and is not discussed fur- 
ther. For simplicity a circular SS orbit is assumed in this report. 


Basic Equations Governing the Problem 

Equations governing the flight of a rocket in an inverse -square gravitational field 
are given by 


G Y e M) ~ 

• u r + e f 

(3) 

r 3 m 

r = V 

(4) 

m = -0(0 

(5) 

i = 1 

(6) 

A A 

f- f - 1 = 0 

(7) 


In these equations, r, V, and m represent instantaneous radius, velocity, and mass, 
respectively; V is the engine exhaust velocity; 0 is the mass flow rate, which may be 
a function of time; and f is the unit thrust direction, which must be determined. The 
superscripts - and * represent vector and unit vector quantities, respectively; and 
' represents a derivative of the particular variable with respect to time. State variable 
C was introduced to remove the explicit time dependence of mass flow rate from the 
equation. This is standard procedure. 

By using a Hamiltonian formulation of the Maximum Principle (ref. 4), the costate 
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equations may be obtained from the so-called Hamiltonian given by 


H = X • V + ju*r+am+r? + y(f*f-l) 


( 8 ) 


In this equation, X, / u, and a are the adjoint variables associated with state variables 
V, r, and m and t is associated with dummy variable The constraint associated 
with the thrust direction is also adjoined to H by y. From equation (8) the time deriva- 
tives of the adjoint variables are given by 

X = -p (9) 

p = iL|~x - 3 (X- r)rl (10) 

3 


a 


VJ - - 
— (X • f) 


m 



where 



m 


o 


(ID 


( 12 ) 


(13) 


The thrust direction that minimizes the Hamiltonian is given by 

A A 

f = X (14) 

and 

y — X (15) 

2m 

State equations (3) to (6) and adjoint equations (9) to (12) make up a set of first-order non- 
linear differential equations which must be numerically integrated in order to obtain a 
solution to the two -point boundary value problem. Optimum thrust direction is along the 
unit X vector given by equation (14). 
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Trajectory Constraints 


In this section, constraints on the problem are discussed in more detail, 
illustrates a simplified trajectory profile departing from a given initial orbit. 


Sketch (b) 
Events a. 



and bj (i = 1 to 3) are also indicated- on the sketch. The outbound burn is initiated at 
time aj = 0 at an optimal point in the initial circular inclined orbit. Thrust is ter- 
minated when the desired mission energy is reached. This is followed by a given coast 
phase until time b 2 . At time (a^ = b 2 ) a retroburn is initiated until time bg is 
reached, followed by the apogee and perigee delta -velocity maneuvers. It may already 
be obvious that a^ represents the time at the beginning of phase i and bj is the time 
at phase i completion. With this convention, constraints governing the problem can be 
written mathematically as follows: 


q! = r(a x ) - rj = 0 

(16a) 

q 2 = V( ai ) - V 1 = 0 

(16b) 

q 3 = rfej) • Vfej) = 0 

(16c) 

A. A 

= h(a^) • l g - cos ij = 0 

(16d) 
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q 5 * ~ h (a ,)- 1 . 

1 2 cos sin ij - o 

% - mfej) - s q 


^7 * r fa 3 ) - F(b 2 ) * o 
q 8 = ^ a 3 > * ^(b 2 ) * 





q Xl = »(b 3 ) - m H e V V 4 V%) 

q 12 = r ft 3 ) • i> 0 * o 

q 13 = V(b 3 ) • h 0 , 0 
q M = b 3 - », + T - t d = „ 
q 15 ' a 3 - bj = 0 
q 16 » ffaj) . Ffbj) , J 

q 17 * V (* 2 ) - V(bj) « ff 

q 18' a 2 + 4t-l, 2 , 0 



= 0 


q 20 = 10 teg) - infbj) » Q 


(16e) 

(I6f) 

( 1 %) 

(m) 

(161) 

(m 

(16k) 

(16*) 

(16m) 

(16n) 

(16o) 

(I6p) 

(16q) 

(16r) 

(16s) 

(m 
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Equations (16a) to (16f) express fixed initial conditions, while (16g) and (16h) express 
the continuity of position and velocity vectors at injection. Equation (16i) gives a fixed 
energy requirement. Equation (16j) is the constraint for fixed declination of the outgoing 
asymptote, (16k) is the constraint on Tug mass required at bg to perform the apogee 
and perigee delta velocity maneuvers and return the Tug empty to the rendezvous orbit. 
The term h Q in equations (16Z ) and (16m) is the unit angular momentum vector of the 
precessed SS orbit to which the Tug must return. These constraints assure that the Tug 
will be in the proper orbit for rendezvous at the end of retroburn. Equation (16n) ex- 
presses the total -trip -time constraint associated with the problem. Equations (16p) to 
(16t) express the constraints at the beginning of payload separation. 

With these constraints, the augmented Hamiltonian to be minimized can be written 
as follows: 


3 20 

H* = m(a 3 ) - m(b 2 ) + ^ H i + Z € i q i 

i=l i=l 

The absolute value of m(ag) - m(b 2 ) represents payload dropped during the coast phase. 
The term Hj represents the Hamiltonians during outbound burn, payload -separation 
coast phase, and retroburn. Constraints q i are adjoined to H* by using the arbitrary 
multipliers e. Auxiliary variational conditions associated with the two-point boundary 
value problem can be derived by using this augmented Hamiltonian. 


Auxiliary Variational Equations 

By using the Hamiltonian Maximum Principle (ref. 5) and the constraints described 
in the previous section, the following auxiliary variational equations must be satisfied to 
maximize payload: 


A 

X(aj) = e 2 V ^ a l ) + e 3 r ^ a l ) + *4^^ + e 5 V r q 5 

(17a) 



fx(a x ) = e 1 r(a 1 ) + tgVfep + e 4 V f q 4 + 

(17b) 

x^) = xfb^ 

(17c) 

p(a 2 ) = fKbj) 

(17d) 


10 



(17e) 


a(a 2 ) = a(bp 


M a g) - X(bg) + €gV(b 2 ) + 


(17f) 


P 

M( a 3 ) = p(b 2 ) + e g — r(b 2 ) + e 10 V-(p 

r 3 


cr(a 3 ) = 1 
cx(b 2 ) = 1 
m(b„) 


(17g) 

(17h) 

(17i) 


X(b 3 ) + [H(bg) - 7<b 3 )] VyT + orCbg) —L Vy AV t + -JL [h Q - h 0 • V(b 3 )V(b 3 )] = 0 

3 


(17J) 


h o -h 0 -r(b 3 )r(b 3 )] =0 

(17k) 

[H(a 2 ) - 7<a 2 )J - [H(b 1 ) - Tfb^J + [h^) - 7(a 3 )] - [H(b 2 ) - -rfbg)] = 0 (111) 

In these equations, represents arbitrary constants. Equations (17c) to (17e) ex- 
press the continuity of adjoint variables at the beginning of the payload -separation coast 
phase, and (17f) and (17g) give the discontinuity in X and p at the end of this coast 
phase. The multiplier a associated with mass is continuous and is equal to 1 at the end 
of the coast phase. Equations (17j) and (17k) must be satisfied at the end of retroburn 
for maximum payload. Gradients of return time T, total delta velocity (AVrp = AV a + 
AVp), and declination cp are derived in appendixes B, C, and D, respectively. Return 
time is defined as the time from retroburn completion to rendezvous with the SS. 


r 1 Ci o 

p( b 3 ) + [H(b 3 ) - T(b 3 )J V— T + a (b 3 ) V F AV t + — — 

e 3 


Evaluation of Discontinuity in State Variables at the End of 
Payload -Separation Coast Phase 

In order to determine the discontinuity in X and p at time a 2 , two arbitrary con- 
stants (e g and e^) must be evaluated in equations (17f) and (17g). First, the arbitrary 
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constant is evaluated as follows: 
Let 


C(t) = A (t) x V(t) + fi(t) x r(t) (18) 

where x represents cross product of two vectors and vector C is a constant of motion 
(C(t) = constant vector). This can be shown by taking the time derivative of C and sub- 
stituting equations (3), (4), (9), and (10) for V, r, X, and Ji, respectively. The time 
derivative is identically zero; therefore, C is a constant of motion. Taking the cross 
product of equation (17f) with V and (17g) with r gives 

C(a 3 ) - C(b 2 ) = e 10 [Vy<p x V(b 2 ) + V-<p x r(b 2 )J (19) 

The constancy of C during retroburn implies that 

C(a 3 )=C(b 3 ) (20) 


Using equations (17j) and (17k) gives 

C(b 3 ) = -h Q x [c 13 V(b 3 ) + e 12 r Ct> 3 >J (21) 

This equation was derived by using the results of reference 6, where it is shown that 
Vyy x V + V— y x r = 0 if y = y £r, V, r • Vj. Since return time and apogee and perigee 
delta velocities are of this form, the corresponding cross products are zero. From 
equation (2 1) it is obvious that 


C(b 3 ) • h Q = 0 (22) 

A, 

Taking the dot product of equation (19) with h Q , substituting equation (22), and solving 
for e 10 give 


CAVA 

10 jvy <p x V(b 2 ) + V-(p x r(b 2 )J • h Q 


(23) 


Given e and X(b 2 ), let 


X - X(b 2 ) + CjQ Vy (p 


(24) 
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With this definition, equation (17f) becomes 


X(a 3 ) = e g V(b 2 ) +X (25) 

In order to determine the discontinuity in X and €q must also be evaluated. 
This may be done as follows : Substitute equations (17c), (17d), (17f), and (17g) into 
equation (17Z), which gives 


0( 3-g) 


V e^ a 3 ) 

m(a 3 ) 


- a(a 3 ) 


- PQoJ 


m(b x ) 


- ofb^ 


= 0 


(26) 


Since /3(a 3 ) = ^(b^ * 0, equation (26) can be divided by /3(a 3 ). Also d=0 during a 
coast phase ( a is constant); therefore, equations (17e), (17h), and (17i) give a(a 3 ) = 
cf(bj) and equation (26) reduces to 


m(a,) 

Mao) = — X(bJ (27) 

J m(b t ) 1 


Substituting equation (27) into (25) and solving for e 9 give 


€g = - 1 ] 

r_ - t 

Jr— - -1 p 

m(a„)X(b 1 ) 

— \ 

2 


/ x- va> 2 r + 


- X* > 

9 V <ty 

l ' 

/ L s J 

m^) 

J 


The choice of + or - sign can be determined as follows: Assume the payload is zero; that 
is, m(a 3 ) = m(b^). Also if the declination is not specified, equation (24) becomes X = 
X(b 2 ). Furthermore, if the coast phase is zero, X(b 2 ) = A(bj). With these assumptions, 
equation (28) becomes 


e 9 = - - 1 — [x(b 2 ) • V(b 2 ) ± X(b 2 ) • V(b 2 )] 

2 

A choice of - sign implies that e 9 = 0 and that the thrust direction is continuous at a 3 , 
essentially pointed along the velocity vector. However, in order to bring the Tug back, 
the velocity must be reduced to below escape velocity, and the thrust direction should be 
opposite to the velocity vector. Therefore the 4 sign should be used in equation (28). 
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Variational Final Conditions at Termination of Retroburn 


In this section, final conditions at the termination of retroburn are derived. Sub- 
stituting equations (17j) and (17k) into equation (8) gives 


[H(b 3 ) - T(b 3 )] 

1- G VyT • r(bg) + VpT • V(bg) 




- a(b 3 ) 


m(b„) 


— — — Vy AV t . r(b 3 ) + v- AV t • V(b 3 ) 

r3 (b 3 i 


However, 


and 


+ ^(bg) 


m(b 3 ) 


- cr (b 3 ) 


G 

r 3 (b 3 ) 


VyT • r(b 3 ) + VjT • V(b 3 ) = jET 


(29) 


G 

r 3 (bg) 


V V 


AV t • r(b 3 ) + V- AV t • V(b 3 ) = 


d AV t 
dt 


These equations are valid for a coast phase (0 = 0). However, the total delta velocity 
required to return the Tug does not change during a coast phase, and d AV^/dt = 0. 
Similarly, the change in return time is equal to the change in time along the coast arc 
but opposite in sign; that is, dT/dt = -1. Substitution of these results into equation (29) 
gives 


a(b 3 ) = 


m(b 3 ) 


(30) 


The arbitrary constants e^g anc * e i 3 * n equations (17j) and (17k) can be evaluated by 
taking the dot product of equation (21) with r (bg) and V(bg): 
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(31) 


and 


€ 13 _ X ^ b 3 ) • h ^ b 3 ) 

V(b 3 )' _ h(b 3 ).h o 


e 12 n( b 3 )-h(b 3 ) 

7^)~~ h(b 3 ).h 0 


(32) 


The remaining two variational conditions can be obtained by taking the dot product 
of equation (17j) with Vy AV T and (17k) with V - AV^,: 


X(b 3 ) - X(b 3 )Vy AV t + |H(b 3 ) - r(b 3 )J VyT 


X(b 3 ) • h(b 3 ) [-* 
" h(b 3 ). h 0 

|p(b 3 ) + x(b 3 )v F av t + [H(b 3 ) - r(b 3 )] v f t 


ii \Uq j r~ A, A, A A — I 

7 r 1 L h o - h o ' v(b 3 )V(b 3 )J U Vy AV T = 0 




p(b 3 ) • h(bg) r A 


h(b~) 


• r a ^ a a n 

[h Q - h Q • r(b 3 )r (b 3 )J V • V F AV t = 0 


3' “o 


(33) 


(34) 


where 


H(b 3 ) - r(b 3 ) = - _iL_ X(b 3 ) • r(b 3 ) + ju(b 3 ) • V(b 3 ) (35) 

r 3 (b 3 ) 

Equations (33) and (34) must be satisfied by the solution of the two-point boundary value 
problem. In order to assure this, they will be used as final conditions to be satisfied 
during the Newton -Raphson iteration. 


Variational Final Condition at Departure from Initial Orbit 

Equations (17a) and (17b) are six equations with five arbitrary constants. This im- 
plies there is a variational equation that must be satisfied. The variational condition 
may be obtained by computing the C vector at a^: 
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C(a.j) — ( e 4^ 3 ~ € 5 ^ 2 ) * b(a^) (36) 

From this equation the variational final condition is obtained by taking the dot product 

A. 

with h(a^): 

C(a 1 ) * h(a^) = 0 (37) 

Initial and Final Conditions 

The initial adjoint variables are unknown and must be guessed in solving the two- 
point boundary value problem. As it is shown in reference 7, these variables can be re- 
placed by physically more meaningful parameters, such as vehicle pitch attitude if/, 
pitch rate if/, yaw attitude 6 , yaw rate 6 , and X(a^) and X(a<) instead of the compo- 
nents of X(a^) and p(a^). With these changes in variables, the free initial conditions 
and associated final conditions are as follows: 


Initial conditions 

Final conditions 


Ma x ) 

r(b 3 ) • h o = 0 

(38a) 

^(a x ) 

V(b 3 ) • h Q = 0 

(38b) 

6 (a 1 ) 

<p[r(b 2 ),V(b 2 )J - = 0 

(38c) 

8 (a x ) 

Equation (33) 

(38d) 

X(a t ) 

cr(a 3 ) = cr(b 2 ) = 1 

(38e) 

X( a l) 

Cfej) • h(a^) = 0 

(38f) 

b l 

E[r(b 2 ),V(b 2 )] -E d = 0 

(38g) 

b 3 

b 3 " a i + T “ t d = 0 

(38h) 

m p = m(a 3 ) - m(b 2 ) 

AV t /V p 

m(bg) - m^e 1 e = 0 

(38i) 

u(a x ) 

Equation (34) 

(38j) 


16 



Condition (38e) can be omitted from the iteration since the entire problem can be scaled 
by X(a^), and this condition is satisfied by scaling the problem. Conditions (38g) and 
(38h) are removed from the iteration by initiating the payload -separation coast phase 
when the desired energy is reached and terminating the retroburn when the desired total 
trip time is satisfied. These two conditions are used as cutoff variables during the it- 
eration. Final condition (38f) is used to compute A(a^), as shown in reference 5. With 
these simplifications, the number of final conditions to be satisfied is reduced to six, 
namely, equations (38a) to (38d), (38i), and (38j). 

In order to solve the two -point boundary value problem, a finite difference Newton- 
Raphson iteration technique is used. The iteration is terminated when the percentage of 
predicted payload change and the absolute value of the normalized desired final condi- 
tions are less than 10“^. 


RESULTS AND DISCUSSION 

The method described in the previous sections was applied to a typical cryogenic 
Tug configuration for a specific set of planetary target conditions. The results of that 
analysis are presented in this section. Performance and weight data for the Tug con- 
figuration are given in table I. Specific impulse has been reduced by 2 percent to account 
for flight performance reserves. Engine specific impulse corresponds to an uprated 
version of the current Centaur RL-10 engine. Combined weight of the Tug and payload 
is assumed to be constant (28 622 kg) and consistent with SS payload weight capability. 
Tug propellant load is varied to maintain this fixed initial weight. The baseline mission 
characteristics are presented in table n. Tug departure orbit is a 18 5 -kilometer - 
altitude circular orbit inclined 28. 5°. The Tug is free to depart from this orbit at any 
point. The departure point is chosen to maximize the injected payload. 

Planetary missions most often are defined by specifying the mission energy and the 
right ascension and declination of the outgoing asymptote. The right ascension for a 
specific mission can be satisfied by selecting the proper SS launch time, as discussed 
earlier. Therefore, this constraint was not included in the analysis, and an optimum 
pseudo right ascension was obtained for each declination. Declination of the outgoing 
asymptote (often referred to as DLA) is parametrically varied from -40° to +40°. A 
mission energy of 12 km /sec was selected as a representative energy typical of an 
inner -planet (Mars and/or Venus) mission. Total Tug trip time for the baseline mission 
is 1 day, and payload-separation coast time is assumed to be 10 minutes. 

Space Shuttle orbital inclination, total trip time, payload -separation coast time, and 
engine thrust are examined later in this section to determine their effects on Tug 
performance. 
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In order to explain the results that follow, it is helpful to recall some well-known 
facts from orbital geometry and previous analyses of planetary missions for expendable 
stages. If the nodal precession of the SS orbit is assumed to be zero and a spherical 
Earth model and planar trajectories are used (no nodal precession during the flight), for 
a fixed initial SS orbit, the same optimum payload can be delivered to all DLA's whose 
absolute value is less than the initial SS orbital inclination. Also there are two distinct 
solutions for each DLA whose absolute value is less than the SS orbital inclination. As 
the DLA approaches SS orbital inclination, the two solutions tend to coincide. At a DLA 
equal to SS orbital inclination, they degenerate into a single solution. DLA’s whose ab- 
solute value is greater than SS orbital inclination can be reached only if yaw maneuvers 
are performed to increase inclination to or beyond the absolute value of the desired DLA. 
These yaw maneuvers reduce performance rapidly. Therefore, for these high DLA’s, 
maximum performance is obtained when the orbital inclination is equal to the absolute 
value of the DLA; and there is a single solution for each DLA. 

When nodal correction associated with the present problem is included, small plane - 
change maneuvers must be performed to correct for SS orbital precession. This will re- 
duce Tug performance, but it should not substantially affect the orbital geometry. It is 
expected that there will be two solutions for DLA's whose absolute value is less than the 
SS orbital inclination and a single solution for DLA’s whose absolute value is greater than 
the SS orbital inclination. 

Another important point to recall is that the rates of change of the orbital inclination 
and the longitude of the line of nodes vary sinusoidally with orbital position for a fixed 
out-of-plane acceleration. The rate of change of the orbital inclination is maximum when 
the out -of -plane acceleration is applied at the line of nodes. The maximum change in the 
longitude of the line of nodes occurs when acceleration is applied at the antinode. This 
implies that, for a fixed amount of propellant expended, the maximum inclination change 
is obtained when the burn is at or near the line of nodes and the maximum change in the 
longitude of the line of nodes is obtained when the burn is added at or near the antinode. 

It is also important to note that an out -of -plane acceleration that tends to increase incli- 
nation near the line of nodes will move the line of nodes retrograde when applied ahead of 
the node and posigrade applied past the line of nodes. Since the line of nodes for the pre- 
sent problem must be moved retrograde approximately 8°, for the DLA's where the out- 
bound burn traverses the line of nodes, nodal shift during this burn is expected to be 
small. 

Results obtained for the baseline Tug and baseline mission are presented in figures 1 
to 6 as a function of DLA. Three SS orbital inclinations were analyzed: 28. 5°, 33. 5°, 
and 38. 5°. The minimum inclination that can be reached with a planar SS launch from the 
Eastern Test Range is 28. 5°. Since the SS does not have substantial dogleg capability, 
lower inclinations were not evaluated. 
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Figure 1 gives optimum pseudo right ascension as a function of DLA. As is ex- 
plained in an earlier section, Assumptions and Simplifications, pseudo right ascension 
permits the user of the data to satisfy specific requirements on the right ascension of the 
outgoing asymptote by selecting SS launch time. By subtracting the value of pseudo right 
ascension for the desired DLA from the required right ascension, the longitude of the as- 
cending node for the SS orbit may be determined. Since the longitude of the ascending 
node for a given SS orbital inclination moves through 360° every day, the needed longi- 
tude of the ascending node may be obtained simply by selecting the correct SS launch 
time. It was mentioned earlier that the SS deploys the Tug and payload into a low Earth 
orbit and that the Tug is free to depart from this orbit at any point. The point of depar- 
ture was determined in solving the two-point boundary value problem; it is presented in 
figure 2 as a function of DLA. The Tug departure argument of latitude is defined as the 
angle measured in the initial SS orbital plane relative to the ascending node of this orbit. 
Positive argument of latitude implies that the Tug departs from the SS orbit past nodal 
passage. The shape of the escape hyperbola is determined by perigee radius, which 
changes little with DLA, and by mission energy, which is constant. Therefore, once a 
trajectory is optimized for a given DLA, the shape of the escape hyperbola does not 
change significantly as DLA is varied. As DLA is changed, the plane containing the es- 
cape hyperbola is essentially rotated about the angular momentum vector of the initial SS 
orbit until the correct geometry is obtained. 

As the value of DLA is varied between the negative and positive values of SS orbital 
inclination, the argument of latitude and pseudo right ascension sweep through nearly a 
180° arc. Recall the earlier discussion on the existence of two distinct families of solu- 
tions for DLA’s whose absolute value is less than SS orbital inclination. Figures 1 and 2 
show only one family. The remaining family of solutions covers the remaining 180° in 
argument of latitude and pseudo right ascension. The existence of these solutions was 
demonstrated numerically by taking the converged initial conditions for a given DLA, 
changing the signs of initial yaw attitude and yaw rate, and adding 180° to the departure 
argument of latitude. With these initial conditions, the trajectory was integrated and 
without any iteration yielded a converged case satisfying all necessary boundary condi- 
tions for a DLA of the same magnitude but opposite sign. The pseudo right ascension for 
these solutions is equal to the pseudo right ascension obtained for the first solution plus 
180°. The nodal correction during outbound burn, orbital inclination at injection, Tug 
payload, and delta -velocity losses associated with the second family of solutions are mir- 
ror images of the first family about a DLA of 0°. Since the second family of solutions 
can be obtained by using the data presented and following the preceding discussion, they 
are not given herein to avoid confusion. 

Figure 3 shows the variation of orbital inclination at payload separation; figure 4 
presents the change in the longitude of the ascending node during outbound burn. In the 
following discussion, figures 2 to 4 and the preceding discussion on changing line of 
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nodes and orbital inclination are used to explain the observed trajectory characteristics 
for the range of DLA’s considered. For high negative DLA’s, the outbound burn starts 
well past the line of nodes, at about 80° argument of latitude (fig. 2). The outbound burn 
arc for the baseline Tug and mission is about 100°. Therefore, the Tug burns through 
maximum latitude, which is an optimal region for changing the longitude of the line of 
nodes. Almost all nodal change is done during this burn (fig. 4). 

Notice that the nodal change during the outbound burn decreases for those negative 
DLA’s whose absolute value is greater than the SS orbital inclination. In order to reach 
these DLA’s, orbital inclination must be increased beyond the initial SS orbital inclina- 
tion (fig. 3). Since the Tug must return to the SS orbit, the same inclination change must 
be made during retroburn, and it becomes more efficient to perform a larger nodal shift 
during this burn. As DLA is increased, the outbound burn moves closer to the line of 
nodes (fig. 2); it becomes less efficient to change the longitude of the line of nodes; and 
nodal change decreases. Also, to move the node retrograde for these cases, the SS 
orbital inclination must be decreased (fig. 3). For DLA's between 0° and 16°, the out- 
bound burn traverses the line of nodes, and nodal change for these cases during outbound 
burn is quite small (fig. 4), as expected. As DLA is further increased, the outbound 
burn begins well ahead of the node and nodal change increases. Also, SS orbital inclina- 
tion increases (fig. 3) and moves the node in a retrograde direction. For DLA's larger 
than the initial SS orbital inclination, the inclination during the outbound burn must be in- 
creased. And again nodal change becomes more efficient during retroburn. The sharp 
corners on the curves for DLA’s nearly equal to the SS orbital inclination in figures 3 
and 4 correspond to the DLA's where the two solutions degenerate into a single solution. 
The reason these corners do not occur exactly at DLA’s equal to the initial SS orbital in- 
clination is that orbital inclination has been changed during the outbound burn to correct 
for nodal precession. 

In figure 5, Tug payload capability is plotted as a function of DLA. Payload is max- 
imum at a DLA of approximately 8° for the three inclinations, and any DLA between -30° 
and 30° can be reached with less than a 100 -kilogram payload loss. Notice the rapid pay- 
load loss as DLA is increased beyond the initial SS orbital inclination. This loss is 
caused by the large inclination change required during outbound Tug burn to reach a 
given DLA (fig. 3) and by the requirement that the Tug return to the SS orbital inclina- 
tion. The inclination change during retroburn must be equal and opposite to the inclina- 
tion change made during outbound burn. Therefore, for a 1-degree increase in DLA, a 
minimum of 2 degrees inclination change must be made. This, of course, becomes very 
expensive in terms of Tug performance. Similarly, the payload is expected to decrease 
rapidly for high negative DLA’s; however, because of numerical sensitivities encoun- 
tered in the Newton -Raphson iteration technique, solutions for these declinations were 
not obtained. For these high negative DLA’s, there is only a single solution; therefore, 
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the performance should be equal to that obtained for the positive DLA’s. If all the solu- 
tions were plotted, the payload curves would be symmetric about a DLA of 0°. 

Figure 6 gives delta -velocity loss as a function of DLA for the three SS orbital incli- 
nations. Velocity loss is defined as the velocity required above ideal impulsive velocity. 
Ideal impulsive velocity is defined as the velocity needed above circular orbital velocity 
at a 185-kilometer altitude in order to reach a given energy. This loss is the result of 
finite thrust level, thrusting along nonoptimal portions of the trajectory, nodal shift 
change, trip time, and fixed payload-separation coast time. Total velocity loss pre- 
sented is the sum of losses during outbound and inbound portions of the trajectory, and 
it remains relatively constant throughout the range of DLA's. At the limiting DLA's, 
both outbound and inbound velocity losses tend to increase sharply, resulting in a rapid 
decrease in payload capability. In order to show the delta -velocity loss due to nodal cor- 
rection, velocity losses for a 90° initial SS orbital inclination (nodal precession of zero) 
are also presented. These cases are discussed later in this section. Notice that with- 
out nodal correction, the velocity losses remain constant. Velocity loss due to nodal 
shift for DLA's between -28° and 28° for an initial SS orbital inclination of 28. 5° is less 
than 75 meters per second (fig. 6(a)). Similarly, the nodal correction velocity losses for 
initial SS orbital inclinations of 33. 5° and 38. 5° are small (figs. 6(b) and (c), respec- 
tively). From the relatively small nodal correction delta velocity required, the assump- 
tion that all nodal correction will be done during the outbound burn and the retroburn 
seems to be a good one. Any increase in payload resulting from introducing an additional 
nodal correction burn would be small. Also, engine startup and shutdown losses asso- 
ciated with a very short burn may substantially reduce or completely eliminate the pay- 
load gain associated with the added nodal correction burn. 

As mentioned earlier , Tug payload capability (fig. 5) is symmetric about a DLA of 
0°. Based on this symmetry, maximum Tug payload capability is given in figure 7. The 
result shown is obtained by departing from the initial SS orbit near the ascending node for 
positive DLA's and near the descending node for negative DLA's. Maximum Tug payload 
capability decreases as initial orbital inclination is increased even though total nodal 
shift decreases. 

In figure 8, Tug payload capability is given as a function of initial SS orbital inclina- 
tion for DLA's of 0°, 10°, 20°, and 30°. Payload capability decreases as initial orbital 
inclination is increased to approximately 45°, then increases with increasing inclination. 
Maximum payload capability occurs at an initial SS orbital inclination of 90°. At this 
point, nodal shift is zero and the. trajectory solution becomes planar (no change in incli- 
nation). Also, payload and velocity losses remain constant as DLA varies for this 
orbital inclination. Based on Tug performance and mission requirement, the best SS 
orbital inclination would be 90°. However, SS payload capability falls off rapidly as 
orbital inclination is increased to 90°. Also, the allowed launch azimuth range from the 
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Eastern Test Range limits SS orbital inclinations to less than 56°. With this SS con- 
straint in mind, for low-DLA missions (DLA range, -30° to +30°), maximum payload 
can be delivered by departing from a 28. 5° inclined orbit. Therefore, this initial orbital 
inclination was selected to investigate the effects of payload-separation coast time, total 
trip time, and engine thrust on payload capability. 

The effect of payload -separation coast time on payload is shown in figure 9. As ex- 
pected, payload decreases as the coast time increases. It is well known that energy- 
change maneuvers can be accomplished most efficiently during high-velocity portions of 
orbits. During payload -separation coast, as radius increases velocity decreases and the 
retroburn becomes less efficient. Therefore, the coast phase should be kept as short as 
possible. The duration of this coast phase (assumed to be 10 min for the baseline case) 
will be determined by the time required to physically separate the spacecraft and turn 
the Tug around for proper orientation to perform the retroburn. Engine exhaust im- 
pingement on the spacecraft is also a concern during retroburn since the Tug is thrust- 
ing essentially opposite to the velocity vector and toward the spacecraft. Thus, suffi- 
cient coast time must be allowed for a safe separation distance to be achieved. 

Effects of total trip time on Tug payload capability are presented in figure 10. The 
total -trip-time constraint is satisfied by varying the retroburn. As total trip time is 
allowed to increase, the amount of energy that must be removed to return the Tug de- 
creases and the payload capability increases. However, the amount of nodal shift in- 
creases as total trip time is increased (eq. (1)), which results in a reduction in payload. 
Therefore, there exists an optimum trip time that gives a maximum payload. Optimum 
trip time for the mission considered varies between 21 and 26 hours in the declination 
range of -30° to 30°. 

Figure 11 shows the effect of Tug thrust level on payload capability for a range of 
declinations. As is expected, the payload asymptotically increases to some maximum 
value as initial thrust increases. These curves were derived by assuming that engine 
weight remains constant as its thrust level varies. In reality as thrust increases, en- 
gine weight increases; and as thrust decreases, so does engine weight. The baseline 
Tug configuration has an initial thrust of 66 723 newtons. For higher initial thrust, en- 
gine weight and consequently Tug hardware weight should increase, decreasing the indi- 
cated payload capability. Based on similar arguments for lower thrust, the indicated 
payload capability should increase. Therefore, there would be a thrust level that gives 
maximum payload. 


SUMMARY OF RESULTS 

In this report, mathematical equations were derived to obtain the maximum payload 
capability of a reusable Tug for planetary missions, taking into account the nodal shift of 
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the initial Space Shuttle (SS) orbit. The analysis assumes that all nodal corrections will 
be done during the outbound burn and the retroburn. The sample case analyzed is an 
inner -planet mission (mission energy, 12 km /see ). The maximum performance of a 
typical cryogenic Tug is shown to be 4790 kilograms for declinations of the outgoing 
asymptote (DLA) between -30° and 30°. It is obtained by departing from a 28. 5° SS 
orbital inclination. Performance varies with DLA; the variation, however, is less than 
100 kilograms for the given range of DLA’s. Declinations outside the preceding range 
can be obtained by injecting the Tug into higher inclination orbits with small performance 
loss. 

The velocity loss due to nodal shift is small for departure from a 28. 5° inclined 
orbit, less than 75 meters per second for declinations between -30° and 30°. The cor- 
responding payload loss is a maximum of about 170 kilograms. 

The total trip time, the elapsed time from Tug departure to return to the SS, that 
gives maximum payload varies from 21 to 26 hours as DLA increases from 0° to 30°. 

Payload-separation coast time was varied from 2 to 20 minutes. Payload decreases 
almost linearly as coast time increases; therefore, to obtain maximum payload, the 
separation coast time should be kept as short as possible. 

Payload is shown to increase with increasing Tug thrust for a fixed engine weight. 
However, engine weight increases with increasing thrust, which would reduce the pay- 
load gain obtained. Consequently, in the real case, there would be an optimal thrust 
level. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, Oct. 16, 1975, 

491-62. 
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APPENDIX A 


SYMBOLS 


a l. time at beginning of phase i, sec 

b. time at end of phase i, sec 

C defined by eq. (18), kg -sec 

DLA declination of outgoing asymptote, deg 

E energy, m /sec 

e eccentricity 

a. 

f unit thrust direction 

2 2 

G gravitational constant of Earth, m /sec 

H Hamiltonian, kg 

H* augmented Hamiltonian, kg 

2 2 

h angular momentum, m /sec 

A 

h Q angular momentum vector of return orbit 
i inclination, rad 

-3 

J oblateness parameter (1. 624x10 ) 

K kappa function, defined by eq. (13), sec 

Z right-handed Cartesian coordinate system components (Z ^ points to initial SS 
orbital ascending node and Zg to North Pole) 

m mass, kg 

P orbital period, sec 

p semilatus rectum, m 

q constraint 

Rg radius of Earth (6. 37816x10^ m) 

r radius, m 

sign sign function 

T return time, sec 

T d trip time, sec 

t time, sec 
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t time from perigee , sec 

At payload -separation coast time, sec 

u argument of latitude, rad 

V velocity, m/sec 

<•» 

Voo direction of outgoing asymptote 

AV apogee delta velocity, m/sec 

a 

AV perigee delta velocity, m/sec 

AV t total delta velocity (AV t = AV a + AV p ), m/sec 

X defined by eq. (24), kg-sec /m 

y dummy variable 

a angle defined by eq. (D2c), rad 

/3 mass flow rate, kg/sec 

y adjoint multiplier associated with thrust direction, kg 
6 yaw attitude, rad 

e arbitrary multiplier 

£ dummy variable introduced to make Hamiltonian independent of time, sec 
77 true anomaly, rad 

6 defined by eq. (D2b) 

\ adjoint multiplier associated with velocity, kg-sec /m 
/i adjoint variable associated with position, kg-sec/m 
a adjoint multiplier associated with mass, sec 

r adjoint multiplier associated with variable kg 

(p declination, rad 

\p pitch attitude, rad 

n longitude of line of nodes, rad 

oj argument of pericenter, rad 

Subscripts: 
a apogee 

D desired 

e exhaust 
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H hardware 
i index 
p perigee 
Superscripts: 

— vector 
A unit vector 
• time derivative, (d)/dt 
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APPENDIX B 


GRADIENTS OF RETURN TIME 

The return time T is defined as the time from the end of retroburn to rendezvous 
with the Space Shuttle (SS). It is expressed mathematically as 

T = I (p 1 + P 2 )-tp (Bl) 

Ct 

where P^ is the period of the intermediate orbit established at the end of retroburn, 
tp is the time from pericenter passage at the end of retroburn, and Pg is the period 
of the second intermediate orbit from apogee to return orbit. The periods can be writ- 
ten as 


/r + r 

P 1 = JL_ (_g P] (B2a) 

1 g^A 2 / 

/r + r 

P = _ 2 J— ( J: EiP ) (B2b) 

2 g»/A 2 / 


where r a and r p are the instantaneous apogee and perigee radii, respectively; and 
r p Q is the fixed perigee radius of the return orbit. 

’ The gradients of T with respect to vector y (where y is either r or V) is given 
by 
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V-r„ = 
r a 


¥-t 


Vp 


hxr.r 2 p V 


V xh _Pp 


V-r n = 
r p 


G Vr 


The partial derivative of time from the pericenter with respect to true anomaly St^/dr) 
is given by 


07j h 


The gradient of true anomaly can be obtained from 


tan„ = h <L v >. 


h 2 - Gr 


The gradients are 


_ „ _ h cos 2 T) - h 2 + Gr 

V - -t r - 

h* - Gr h - Gr 




^ - h cos 2 v Tr h 2 + Gr / r • V \ 
h 2 - Gr h 2 - Gr ' h ' 


V-h + ^V); 
h 2 - Gr 


The gradients of h are given by 


Vyh = h x r 
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and 


Vph = V x h 

From these equations the gradients of T with respect to radius and velocity can be 
computed. 


(B8b) 
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APPENDIX C 


GRADIENTS OF TOTAL DELTA -VELOCITY CHANGE 

The total impulse delta-velocity change is the sum of the absolute values of apogee 
and perigee delta -velocity changes; that is, 

AV t = |AVj + |AV p | (Cla) 


where 


2Gr 


AV a = 


P,D 


2Gr, 


f r a^ r a + r p,D ) V r a (r a +r p ) 


(C lb) 


AV =■ 
P 


2Gr 


2Gr 


a,D 


(Clc) 


r p,D^ r a + r p, D^ Y^.D^D^p,^ 

Again let y stand for either radius or velocity; then the gradient of AV T is given by 

V- AV t = sign(l , AV a )V- AV a + sign(l, AV p ) V- AV p (C2) 

The gradient of apogee delta velocity is given by 


V- AV a = 




afra-VD^p.P^a + VD 1 ., 

p,P r|(r a + >- piD ) ? 


2Gr 


f a< r a + r p> 

1 

r P (2r a + V 

i 2Gr P 

_( r a + r p) 

o y p 2 

i 2 4 

( r * * r pf 


V-r 


y a 


(C3) 


Similarly, the gradient of perigee delta velocity is given by 
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1 


(C4) 


V- AV = 

y p 


.. J r p,D (r a+ r p,D> 
V 2Gr 


( r a + r P.D) : 


V-r 

2 y a 


Gradients of apogee and perigee radii with respect to radius and velocity are given by 
equations (B4). 
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APPENDIX D 


GRADIENTS OF DECLINATION 

Declination is defined as the angle between the outgoing asymptote and its projection 

A A, 

in the equatorial plane. Let V M be the direction of the outgoing asymptote and l^ = 

(0, 0, 1) be a constant vector pointing to the North Pole. Then the declination is given 
by ' ” 

sin <p = Z 3 • V M (Dl) 

The outgoing asymptote can be expressed as 

A A A 

= cos a r + sin at e (D2a) 

where 

0 = h x r (D2b) 

and 

a = cos” 1 - — - 77 (D2c) 

e 

where 7 / is the true anomaly and e the eccentricity at injection. Let y represent 
either radius or velocity; then 

V-<p = £2£_^ [(-tan a r + £>) • L V-a + V-(L • r) + tan a V-(L • e)l (D3) 

y COS (p *- 0 y y 0 y 0 j 

The gradient of a is given by 


V ’ -V — 7= V 

e je^ - 1 

Gradients of true anomaly can be computed from equations (B7) and (B8). The gradients 
of eccentricity are given by 


V-e = JL (2EV_h + hV-E) 

y o y y 

eG z 


(D5a) 
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The gradients of h are given by equation (B8). The gradients of energy are as follows: 



l> 

II 

w 

(D5b) 


V-E = — r 

r 2 
r 

(D5c) 

In order to compute the gradients of declination, it remains to evaluate the gradients of 

• r and l 3 ■ 6. 

This can be done from the following equations: 



Vy (l 3 * r ) = 0 

(D6a) 


a a < Pa a a a "1 

V ! 3- r)= ;l i 3-('3- r)r J 

(D 6 b) 


A A 

A A 1“ A A A A^ l O * 0 

V i 3'«)^l'3-«3- r)r ] — y-V 

(D 6 c) 

v F<V «>- 

n \T A A 1 A A [ A XT 1 A A A A A A ■* A A 

61 ( h .V) -±(Zg- e) r I (Z 3 - r)V + (r- V)Z 3 +I (Z g . eJVph 

(D6d) 
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TABLE I. - CHARACTERISTICS OF 


BASELINE TUG 


Initial weight (including payload), kg . . . 28 622 


Burnout weight, kg 2767 

Engine thrust, N 66 723 

Engine specific impulse, I , sec 460 

Performance reserve, percent of I 2 


TABLE n. - CHARACTERISTICS OF 
BASELINE MISSION 


2 2 

Vis-viva energy, km /sec 12 

Range of declinations, deg -40 to +40 

Total Tug trip time, days 1 

Payload -separation coast phase, min ... 10 

Initial orbital altitude, km 185 



Figure 1 - Pseudo right ascension of outgoing asymptote measured in equatorial plane 
counterclockwise from ascending node for baseline mission and baseline Tug. 
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Figure 11. - Tug payload capability as function of Tug engine thrust for constant engine 
weight, 28.5° Space Shuttle orbital inclination, baseline mission, and baseline Tug. 


NASA- Langley, 1976 E — 8498 


41 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON. D C 20546 


OFFICIAL BUSINESS 

penalty FOR PRIVATE USE $300 SPECIAL FOURTH-CLASS RATE 

BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 

4SI 



POSTMASTER 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute . to the expansion of human knowl- 
edge of phenomena m the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof" 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge 

TECHNICAL NOTES Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge 

TECHNICAL MEMORANDUMS 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge 


TECHNICAL TRANSLATIONS Information 
published in a foreign language considered 
to merit NASA distribution in English 

SPECIAL PUBLICATIONS Information 
derived from or of value to NASA activities 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies 

TECHNOLOGY UTILIZATION 
PUBLICATIONS Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 




